Defects in nematic membranes can buckle into pseudospheres 
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A nematic membrane is a sheet with embedded orientational order, which can occur in biological 
cells, liquid crystal films, manufactured materials, and other soft matter systems. By formulating 
the free energy of nematic films using tensor contractions from differential geometry, we elucidate the 
elastic terms allowed by symmetry, and indicate differences from hexatic membranes. We find that 
topological defects in the orientation field can cause the membrane to buckle over a size set by the 
competition between surface tension and in-plane elasticity. In the absence of bending rigidity the 
resulting shape is universal, known as a parabolic pseudosphere or a revolved tractrix. Bending costs 
oppose such buckling and modify the shape in a predictable manner. In particular, the anisotropic 
rigidities of nematic membranes lead to different shapes for aster and vortex defects, in principle 
enabling measurement of couplings specific to nematic membranes. 



I. INTRODUCTION 

By "nematic membrane" we refer to any flexible sheet 
incorporating ordered rod-like constituents. For exam- 
ple, thin films of smectic-C liquid crystals are nematic 
membranes [TJ[H|3J|1]. Also, recently developed sheets 
of carbon nanotubes have nematic character [SJ |5] . Ne- 
matic order arises in lipid membranes with inclusions [7] 
and in the cell cytoskeleton, e.g. during mitosis [8j. Inter- 
estingly, in vitro experiments on mixtures of cytoskeletal 
filaments and protein motors observe topological defects 
(asters and vortices), which spontaneously self-organize 
into a variety of patterns [9l[10l[TT]. These experiments, 
and related simulations, use flat geometries with various 
boundary conditions ^12, il3j . Similar topological defects 
influence the shapes of real cells. For example, cells of the 
alga Bryopsis sprout branches out of vortex-shaped de- 
fects that appear naturally in their cell wall of cellulose 
microflbrils |14j . To take a step toward understanding 
such living and in vitro systems, we consider equilibrium 
shapes around defects in deformable nematic membranes. 

We show that topological defects can buckle the mem- 
brane. This has similarities to two other systems. One is 
bulk nematic liquid crystals, which buckle into the third 
dimension around defect lines |151 1161 117] in a manner di- 
rectly analogous to the shapes we flnd. A second example 
is provided by deformable triangular latices, which have 
been studied extensively in the theory of two-dimensional 
melting. While the physical picture is different, the 
model energy is equivalent to a nematic membrane with 
isotropic elastic constants. Disclination defects culminat- 
ing in a site with five or seven bonds (instead of the usual 
six) can lower their energy by buckling [151 IT^ EUl ET] . 
When draped over curved surfaces, collections of such de- 
fects arrange in specific patterns [22j [23j [24l [25] . If sur- 
face tension is neglected, five-fold disclinations assume 
an approximately cone-shaped form [19 . 

In contrast to the above cases, competition between the 



cost of surface area and rod misalignment determines the 
shape of the defects we consider. When bending rigid- 
ity is neglected, we find that topological defects deform 
membranes into a simple universal shape known as a 
parabolic pseudosphere [23^, ^7\. The size (height and 
extent) of this universal form is governed by the ratio 
of surface tension to in-plane elasticity. The inclusion of 
bending rigidity opposes this puckering. If the bending 
cost is small, the singularities at the tip and rim of the 
defect become smoother. The logarithmically diverging 
tip of the parabolic pseudosphere is replaced by the finite 
height of an elliptic pseudosphere [27j, and the sharp rim 
is replaced by an exponential falloff with a length scale 
related to rigidity. Higher bending costs completely elim- 
inate the buckling instability. The anisotropic elasticity 
of nematics singles out specific defect orientations (asters 
and vortices) ; and corresponding anisotropics in bending 
rigidity lead to different length scales for their shapes. 

The rest of the manuscript is organized as follows: In 
Sec. [n] we describe the free energy of a nematic mem- 
brane using tensor contractions from differential geome- 
try. This provides a compact formulation applicable to all 



deformations, including highly curved shapes. In Sec. Ill 



we describe vortex and aster defects, derive shape equa- 
tions for radially symmetric configurations, and solve 
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them to find the buckled defect shapes. Section IV pro- 
vides a summary and indications for future research. In 
Appendix [X] we study filament orientations in fixed ge- 
ometries, which may provide other ways of measuring the 
nematic membrane parameters. In AppendixjB] we check 
the linear stability of the buckled defect shapes. 



II. ELASTIC FREE ENERGY OF NEMATIC 
MEMBRANES 

Using differential geometry to describe a two-surface 
in three-space, we construct a power series expansion for 
the free energy by selecting a linearly independent set of 
scalar contractions of the surface tensors. For a surface 
described by an embedding vector X(a^,a^), one con- 
structs tangent vectors, ti = diX, by taking derivatives 
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of the embedding vector with respect to its two parame- 
ters. The metric tensor is then gij = ti ■ tj . The covariant 
derivative is defined such that D^gj^ = 0. The curvature 
tensor is constructed from covariant derivatives of the 
tangent vectors as Kij — {Ditj) ■ N, where iV is a surface 
normal. One must choose a side to define the sign of N. 

In the principle directions basis, Kat = ^ ^^q^^ l/^R ^ ' 
where Ri are the radii of curvature [28\ . 

A unit-magnitude tangent vector field T = T^U repre- 
sents the nematic particles. At constant filament density, 
the magnitude T'T; = 1 is fixed and only its orientation 
changes |55j . Nematic symmetry implies invariance un- 
der T' — > — T* [SB]. A complete set of scalars up to second 
order in derivatives is 

+ ^{TlK,,Ti-H^)' 

+ ^(T^K,,Ti-H^y . (1) 

This free energy density must be integrated with a sur- 
face area element dA = ^/g (Pa, where g is the deter- 
minant of the metric. The weighted antisymmetric ten- 
sor jij — eij rotates one-tensors by 7r/2, such that 
T\_ = [28j. Each term is manifestly positive, so 

stability demands that the moduli be positive. In the 
remainder, we consider reflection symmetric, non-chiral 
membranes without spontaneous curvatures i?||,_L,x- 

Unlike parameterizations used to study nematic mem- 
branes near the hexatic fixed line [IHl ISH], this set of 
scalars cleanly delineates the anisotropic bending ener- 
gies that make nematic membranes unique. Creating 
more surface area costs a [3TJ [32]. In-plane splay and 
bend cost Ki and K^, which are the two-dimensional 
analogs of the bulk nematic Frank constants [H 155] . 
Membrane curvature in the direction of the local fil- 
ament orientation costs Curvature perpendicular 
to the filaments costs k±. These out-of-plane bending 
terms are the anisotropic analogs of the the Canham|34j- 
Helfrich[35] bending rigidity. Saddle curves cannot be 
constructed from the other two out-of-plane bending 
terms and incur an independent energy cost of Kx. The 
square of the chiral scalar TKT± [29j [36] is non-chiral. 
The underlying membrane has a fiuid character in that 
the particles can rearrange in the surface without stretch- 
ing or shearing costs. 

Compared to the splay, bend, and twist of bulk ne- 
matics, nematic membranes have additional freedom that 
comes from relaxing a constraint: instead of three fields 
constrained to a unit vector, the nematic membrane con- 
strains only two fields to a unit vector and allows a third 
field to range freely in describing the membrane's local 
deviation from flatness [57] . 



In a system of motor proteins pulling on cytoskeletal 
fllaments, Ki would be proportional to motor density, 
which we assume to be uniform, and k± would be deter- 
mined primarily by the bare membrane's isotropic rigid- 
ity. Filament rigidity would influence both Kt, and 
See Appendix [A| for comments on Kx. 

Perturbative RG near the hexatic fixed line f3D] 
shows that thermal fluctuations reduce weak anisotropy, 
i.e., the three quantities Kx — K|| — K|| ~ and 
Ki — K3 fade at long distances, so that only the hexatic 
membrane energy remains, and 

^he.aUc = CT+ ^{D,T,){D'T^) + {K ~ H^f , (2) 

where Ka — \{Ki + K^) and k = + Under 
further rescaling, k ^ Q and Ka is unrenormalized. Note 
that while the hexatic energy takes its name from the six- 
fold symmetry of triangular lattices, any n-atic symme- 
try with ri > 3 restricts Ki — K3 and k± — — 2kx- 
For polar (n — I) or nematic (n — 2) membranes, the 
isotropic approximation is an important limiting case at 
one extreme of a phase diagram that deserves further 
attention. 

Estimates of the thermal persistence length, ^t, of 
weakly anisotropic rigid membranes indicate an exponen- 
tial form log oc K [SUJ [371 [M] • Modest changes in k can 
thus sweep the persistence length from small values up 
to thousands of times the short-distance cut-off [39] • Ef- 
fects unique to the nematic membrane can then appear in 
patches of material smaller than this persistence length. 



III. BUCKLED DEFECTS SHAPES 

In the nematic phase, the rod orientation varies slowly 
throughout most of the material. However, at particular 
defect points, the orientation may be undeflned, because 
rods at neighboring locations point in all directions. The 
topological charge of a defect is the number of times that 
the orientation rotates through 27r as the coordinate an- 
gle 9 sweeps through 27r. Different patterns appear for 
integer, half integer, and positive and negative charges. 
The defect depicted in Fig. [T] is radially symmetric, and 
is rotated by a uniform angle ^ with respect to the radial 
vectors. 

In the limit of isotropic rigidity, +1 defects with any 
radially uniform ^ have the same energy. The symme- 
try is removed by the anisotropic moduli in a nematic 
membrane, which distinguish asters = 0) and vortices 

= '''/2)- The energy of such a planar defect as a func- 
tion of ^ is 

EplanariO = ^ (^1 COS^ (C) + K3 siu^O) In ^ + i^c(0,(3) 

where R is the size of the membrane, a is a short distance 
cutoff, and Ec is a core energy reflecting the defect's mi- 
croscopic situation inside of the core radius a. 
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FIG. 1: Rod orientations around a general uniform +1 topo- 
logical defect. ^ = corresponds to an aster, and ^ = 7r/2 to 
a vortex. 

For > Ki, asters have lower energy than vor- 
tices and are stable against in-plane deformations. If 
Ki > 2k±, the defect energy is further reduced by buck- 
ling out of flatness to align the filaments in the third di- 
mension. (See Appendix [B] for linear stability analysis.) 
Buckling comes at the expense of creating more area, so 
surface tension sets the size of the deformation. Analo- 
gously, when Ki > > 2k||, vortices are stable and can 
reduce their core energy by tilting the surface axound the 



defect. 

To study this buckling, we minimize the nematic mem- 
brane energy, Eq. (jlj, around fixed aster and vortex ar- 
rangements. For a radially symmetric surface with no 
overhangs, we use the polar Monge tangent representa- 
tion with embedding vector X{r,9). The height above 
the Monge plane is found by integrating the tangent an- 
gle xif) from a base value, so that 

(r cos{9) \ 
r sin(6') . (4) 

This yields a metric with no derivatives and thus lower 
order shape equations. To handle shapes with overhangs, 
such as prolate vesicles [ID], one can parameterize the 
shape by contour length instead of Monge radius. 

The unit vector constraint is enforced by defining the 
angle ^ such that 

i _ ( cos(C) cos(x) \ 
^ - [ sin(e)/r ) ■ 

With this parameterization, the nematic membrane free 
energy becomes 



F„. 



2tt / dr 



V 



cos(x) 



K||(sin2(^) tan^ix) + cos2(^)x'^)+ 

K_L (cos^ (^) tan^ (x) -I- sin^ (^) x'^)+ 

Kx sin^(^) cos2(^)(tan(x) - rx')'^+ 

i^i(cos(0 - rsin(Of')' + -^3(sin(e) + rcos{OC? 



cos(x) 
2r 



(6) 



where Kx = Kx — — k±. We could have written the energy directly in terms of Six by switching from a (TKT)^ 
to a TKKT parameterization as permitted by the unit- vector constraint. Fixed aster or vortex configurations carry 
no energy cost from the term proportional to Kx. Setting to zero the functional derivative of Fnematic with respect 
to X yields a shape equation, which for an aster = 0) becomes 







/.M^ + ^(_2cos(x)x' + rsin(x)x'' - 2r cos(x)x") + ^ f 1 
cos-^^X) 2 2 \ 



1 \ sin(x) Ki sin(x) 



cos2(x) / r 



(7) 



For fixed vortices, the same shape equation holds after 
switching the coefficients k_l K|| and Ki ^ K^. 

For any membrane (hexatic or nematic) without stiff- 
ness {{k} = 0), defects have a simple universal shape 
resulting from the competition between the in-plane mis- 
alignment cost and surface tension. The misaligned rods 
near the defect core can align by bending out of the plane 
into the third dimension, at the cost of increasing surface 
area. The optimum tangent angle is given by the simple 
formula 

, , 2ar^ r 



where ro = \jK\l2a is the distance outside of which sur- 
face tension dominates and flattens the surface. Integrat- 
ing the angle gives the universal shape 




which approaches vertical at r = where the height is 
logarithmically divergent. This may be regulated by a 
cut-off, such as the membrane thickness. As a reference, 
at half the rim radius: /i(ro/2) « 0.45 rg. In a hexatic 
membrane K\ = K-j, , so asters and vortices have the same 
radius. In a nematic membrane, asters and vortices have 
different radii; the lower energy defect also has smaller 
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size. 

This shape, Eq. ([9]), is known as a parabolic pseudo- 
sphere or antisphere, because it has constant negative 
Gaussian curvature equal to —l/r^ [57]. It is also 
known as a tractrisoid, because it is half the surface of 
revolution generated by revolving a tractrix about its 
asymptote [H] . The tractrix is the path of an object be- 
ing dragged by a string of constant length along a straight 
line that does not intersect the object. Leibniz likened 
this problem to a dog owner dragging his hound by its 
leash and named the solution hundskurve. The hund- 
skurve has been studied by Huygens and others . This 
construction makes it clear that the distance to the axis 
along the line tangent to any point on the surface is con- 
stant, i.e. the leash length is tq. These shapes of constant 
negative curvature are also known in quantum gravity as 
solutions to classical Liouville theory [43] . 

This simple shape has singularities at the origin and 
at the rim r = ro, which are modified by the mem- 
brane bending rigidities, k± and respectively. Setting 
K| I = removes all derivatives of x from the shape equa- 
tion, so a simple rearrangement provides the solution, 



cos(x) 



' 2ar^ + Kj_ 
Ki - K± 



\ 



1 



(10) 



where ri = y/ {K 

C = K±/{Ki 



1 — 2Kj_)/2cr is the new rim radius, and 
2k±) is related to the now finite slope at 
the tip. For sufficiently large Ki, the surface puckers out 
of the plane for r < ri , with a profile 



h{r) 



tan 



zbcos^i 



Kl - K_L 



dr(ll) 



(12) 



For Kl < 2k± or for ri < r, this solution is not real, so 
X — becomes the only solution to the shape equation. 

Equation ( 12 ) is a complete elliptic integral of the sec- 
ond kind 44j~^We change variables u y/csmh{u) to 
obtain 



h{r) = ±riVl 




ccosh^(u) 



du , (13) 



where the integration ranges from sinh~^(r/(ri-y/c)) to 
sinh~^(l/i/c). In the study of surfaces with constant 
Gaussian curvature, Eq. ([l3| is a familiar expression for 
an elliptic pseudosphere '27J. Figure |2] shows an example 
elliptic pseudosphere. 

The bending rigidity cuts off the logarithmically di- 
verging tip near the core. Near the origin, the elliptic 
pseudosphere is approximately cone-shaped with slope 



h{r) 



±- 



(14) 



If the term in Eq. ([To| were not present, the shape 
would be a cone. Unlike the cone, pseudospheres have 
constant Gaussian curvature. 



1 



R1R2 



det 



1 



det 



N ■ dd^ 



rlil + c) ' 



(15) 
(16) 



where N is the unit vector proportional to ti x t2 and we 
have carried through the computation after inserting the 
coordinate tangent vectors for the elliptic pseudoshere, 



/ 



ti 



cos{9) , sm{6) , 



V 



1 - 


( r\ 


1 


c + 




2 

1 



t2 = (-rsin(6l),rcos(6'),0) 



(17) 
(18) 




FIG. 2: An example of buckled shape for ri — 1 and c = 1. 

Even with finite kj_ the above shape retains a cusp- 
like singularity at the origin. We may well question how 
the singularity is modified by inclusion of cut-offs and 
higher order terms. A simple short-distance cut-off, a, 
can be introduced as the radius of a hemispherical or 
similar cap over the singular point at the origin. The 
curvature energy density ^ 1 /a^ integrated over the cap's 
area ~ leads to a finite energy. We can then regard 
this as a benign singularity that adds a constant to the 
defect core energy Ec in Eq. 

Substituting Eq. (10 1 into the full shape equation, 
Eq. (j7|, leaves a term proportional to both and to r, 
so the elliptic pseudosphere is expected to remain valid as 
r — > near the core. The situation at the rim is very dif- 
ferent: Designating the distance from ri by e = 1 — r/ri, 
one sees that xoc-^/ease— >0^ and is zero immediately 
outside this radius. The abrupt rim would cause the en- 
ergy proportional to K|| to diverge, so when K|| > 0, the 
defect shape must be different. Since x tends to zero 
away from the core, we linearize the shape equation for 
small X E^nd x' to 

~ (^(ifi-2^^)i-2ar^x + 2/«||(x' + r-x") -(19) 
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After changing variables to e and redefining x — > x (e) to 
be a function of e, the hnearized shape equation is 



= e(e-2)x- 



2^||(e-l) 
Ki - 2kj_ 



(x' + (e-l)x") • (20) 



Note that the approximation is made for small x ^^'^ 
x' , and e need not be small. For real- valued x(e)i this 
equation is solved by modified Bessel functions of the sec- 
ond kind with imaginary order. The order and argument 
both diverge with vanishing as 



X(e) cx K 



ri 



with 



iKi - 2k_l 



2k 



II 



(21) 



(22) 



This solution decays exponentially and has no zeros for 
ri < r. Since our parameterization does not handle over- 
hanging surfaces, x is limited to the range (— 7r/2, 7r/2). 
Thus, for a given value of v, the amplitude must be such 
that the solution stays in this range. For of order one 
and larger, an amplitude of unity yields a x that is suffi- 
ciently small for ri < r that the linearized shape equation 
is valid. It approaches zero asymptotically, so the rim 
radius at which x = shifts to infinity. The asymptotic 
form of Eq. (pTl is gSlllS] 



for (j/, r) y oo 



(23) 



which shows that bending rigidity introduces a new 
length scale 



r2 



2ri 
3 ly 



2^11 
3cr 



(24) 



As shown in Fig.|3] this solution for x oscillates sharply 
near the core, which invalidates the small x' approxi- 
mation. In this region, nonlinearities take over and the 
shape crosses over to the elliptic pseudosphere. As K|| ap- 
proaches zero, the Bessel function becomes flat for ri < r 
and oscillates rapidly inside this radius, thus restoring the 
abrupt rim. 



IV. CONCLUSIONS AND OUTLOOK 

We predict that topological defects can buckle ne- 
matic membranes into elliptic pseudospheres with ex- 
ponentially decaying rims. Let us explore this possi- 
bility in the case of a freely suspended thin film of 
smectic-C liquid crystal. The observation that these 
films maintain orientational order at room temperature 
suggests that Ki (and K^) are at least of the order of 
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FIG. 3; The tangent angle (not the height) as a func- 
tion of radius in units of ri. Solid curve: The Bessel 
function solution for x{^) for i/ = 1. Dashed curve: 

cos-i (^^((r/n)' + c)/(l + c)) 

gest matching in the crossover. 



for c = 2.5 chosen to sug- 



4 X 10 erg. If we assume a surface tension 



typical of bulk liquid crystal interfaces, a 10 ergs/cm^, 
and ignore bending rigidities ({k} — 0), then from 
K ^ fcsTloom, we would estimate a rim radius tq = 
yjKjla K, 1/2 nm, which is smaller than a typical film 
thickness and beyond the limits of this coarse grained 
model. To create larger (observable) defects requires ei- 
ther smaller surface tension or larger in-plane stiffness 
than this initial estimate. Some physical systems may 
allow this. 

For example, studies of thin films of liquid crystal of- 
ten observe a small surface tension, because the chemical 
potential for particles in the film is similar to the chem- 
ical potential in the meniscus surrounding the film [31]. 
This reservoir on the edge of the suspended film allows 
the film to increase its area at a low energetic cost. One 
might control the size of buckled defects by manipulating 
the surface tension via this reservoir. 

In addition to having a surface tension smaller than 
our initial estimate, some materials have observed val- 
ues of {K\ one or two orders of magnitude larger than 
room temperature. For example, scattering studies by 
Spector et al on thin films of smectic-C 80SI found large 
values of {K\ and surface tension made small by the 
meniscus [3]. This particular study used a smectic-C tilt 
angle of 32.2° and found K^jKj, = 4.6(±0.4), k\\IK^ = 
3.4(±0.3), Kx/k|| = 5.5(±3.1), = 75(±24), and 

K|| ~ 10~^^ergs. Unfortunately, the large value of Kn/iCa 
prevents buckling. Since Ki/K^ > 1, flat vortices should 
be stable relative to asters. For other tilt angles or other 
materials, one might hope to find lower values of k^^/K^ 
that allow buckling. 

For sufficiently fioppy films, buckled defects could be 
observable via specular reflection or by interferometry 
techniques used to measure the flatness of mirrors. Addi- 
tionally, islands of smectic-C material may provide means 
of manipulating single defects with laser tweezers, al- 
though coupling between the island's multiple smectic 
layers may introduce additional affects [47 \ . 
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While smectic-C 80SI has sufficiently small surface 
tension and large in-plane stiffness, its bending rigidity 
suppresses buckling. That such shapes have not been 
observed so far in other materials may well be an indica- 
tion of the importance of bending rigidity. Since typical 
lipid membranes have k ~ SfesTroom, this is a severe 
constraint. We note, however, that for stiff rods (nan- 
otubes, cytoskeletal filaments) embedded in membranes, 
the rigidities, Kj| and kj^, and corresponding stiffnesses, 
Ki and K^, may well differ by orders of magnitude. The 
challenge remains to obtain estimates of these parame- 
ters for specific microscopic models, and come up with 
an appropriate system for the study of buckled defects. 
Observations of these shapes in nematic membranes may 
provide estimates of the ratios between «;||, Ki, K^, 
and a. Measurements of Kx are possible via other shapes 
as described in Appendix [A] 

If one could control the anisotropic bending rigidities 
individually, one might be able to sweep a nematic mem- 
brane through a sequence of regimes in which different 
types of defects are stable. For example, for < Ki, if 
one could hold K|| fixed while adjusting k± one might ob- 
serve buckled vortices when 2k± < Ki — K^ + 2«;||, and 
buckled asters when Ki — K3 + 2k| | < 2k± < Ki , and 
buckled vortices again when Ki < 2k±. By increasing 
K|| while keeping k± in any of these regimes, one would 
flatten the preferred shape of vortices. Thus, it is possi- 
ble for vortices and asters to prefer buckled or flat shapes 
independently. 

In focusing on shapes of minimal energy, we have ne- 
glected thermal fluctuations. At long distances, thermal 
fluctuations reduce the differences between aster and vor- 
tex defects in weakly anisotropic membranes |30j . In fu- 
ture work, we would like to explore if this is still the case 
in strongly anisotropic membranes, or if thermal fluctu- 
ations can enhance the anisotropy. 

Unlike bulk nematics, nematic sheets often appear with 
naturally periodic boundaries such as closed vesicles. By 
the Poincare-Brouwer theorem j48] , a genus zero nematic 
vesicle must have topological charge of +2. In fact, de- 
fects can burst the vesicle [49^ 1501 EI]- This resembles 
Bryopsis sprouting branches out of defects in its tethered 
nematic cell wall - a topic to which we hope to return in 
the future. 
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APPENDIX A: NEMATIC MEMBRANES IN 
FIXED GEOMETRIES 

While Kx did not contribute to the shape of -1-1 de- 
fects, it affects other geometries. Following de Gennes' 
molecular field argument [52], we impose the unit vector 
constraint via a Lagrange multiplier A and seek energy 
minimizing filament configurations in fixed geometries. 
In a two-bein basis aligned with the principle directions, 
the curvature tensor is diagonal and T' = e^Ta, where 
is a transformation to local coordinates in which q, , — 



at each point 



In the principle two-bein. 



cos(0 
sin(0 



(Al) 



In the following, there is no summation over a or 
a={a+l) mod (2). The functional derivative of 
Eq. ([TJ in the principle two-bein reads 



Ha ^ ^ = -Xia\a')Ta 



(A2) 



= {k^^CI + K^Cl + SxT|(Ci - C^f] (A3) 

-Kida {D,T,)) - Ksetadb {ecdDcTd) 
= Ta {k\\CI + K^Cl + SxT|(Ci - C^f] (A4) 

-K^da (5e^ - A,) Tie 

-Ksebadb (dc^ - A,) T, , 

where Ca are the principle curvatures and Aa — ei ■ 9a 62 
is the spin connection. The three-vectors ia = e\ti form 
an orthonormal basis in the principle two-bein. To setup 
the molecular field equation, one must carry out the 
derivatives and puUout an overall factor of Ta to obtain 
an expression for A that is a function of index a. One 
obtains an equation for ^ by requiring A to be a scalar, 
i.e., to have the same value for both a = 1 and a — 2. 
Solutions to this equation for ^ extremize the energy. 

Considering first Ki = = A'3, the equation yields a 
simple solution for ^, 



cos2(0 



1 



1 



C1 + C2 



HI 



(A5) 



This is only valid with both components of Ta are non- 
zero, so ^ = and ^ = 7r/2 must also be considered in the 
list of possible ^ values. One must check which candidate 
value for ^ minimizes the energy for particular values of 
{k} and the principle curvatures. In the following, we 
list a few special cases. When k = K|| = kj^ 7^ '«x/2 and 
Cl ^ C2, one has that ^ = 7r/4 minimizes the energy if 



> 



(Cl - C2)' 
4C1C2 



(A6) 



Otherwise, the filaments align with the least curved di- 
rection. 

When Kx = K|| + so that Kx = 0, the orientation 
can be found by minimizing the energy with respect to 



7 



TABLE I: 



Angle 


Stability Criterion 


^ =0 


Ci«:|| > C2K± 


? = V2 


C\n± > C2K11 


cos (t,) - (Ci_C2)(«i+K||) 


ClK± < C2KII, 

CiKii < C2«:± 



tion patterns on curved shapes. Substituting the princi- 



ple curvatures for the elliptic pseudosphere into Eq. ( A5 1 
yields an equation for ^ that is not constant, 



(Bl) 



directly, instead of the molecular field equation. The 
result for Ki — — K3 and Ci < C2 is shown in Table |l] 

On a developable surface, i.e. Ci — and C2 7^ 0, 
when Kx > 2min the stable orientation is 

aligned with the uncurved direction. For smaller values of 
Kx , a special intermediate angle is the global minimum. 



cos2(0 



2k 



II 



2 (k|| + k_l - Kx) 



(A7) 



Note that this only occurs when both parallel and perpen- 
dicular bending are more costly than Kx/2. This might 
result from rods that weaken the sheet or have a specific 
texture on the rod's surface. 

For a developable surface, the spin connection is 
zero, so the covariant derivatives become regular par- 
tial derivatives. Thus, on a cylinder, far from 
boundaries, a constant orientation solves the full 
molecular field equation with the gradient terms in- 
cluded. This could allow experimental measurement of 

(2k|| - Kx)/(k|| + K±- Kx)- 

In more general geometries, in-plane splay and bend 
compete with out-of-plane bending in a non-linear PDE, 
which, in principle, can be numerically integrated to fit 
model parameters to vectorized images of a real nematic 
membrane. Computing model parameters from such im- 
ages in the presence of topological defects requires care. 



APPENDIX B: STABILITY OF BUCKLED 
DEFECT SHAPES 

As discussed in Appendix [K\ the relative strength of 
Kx plays an important role in the stability of orienta- 



and thus not the perfect aster (or vortex) that we as- 
sumed when setting up the shape equation, Eq. ([t]). Since 
Eq. (A5l was derived assuming Ki — — K^, the ques- 



tion remains whether the buckled defect is stable to per- 
turbations away from a perfect aster (or vortex) . 

To check this, we construct linearized evolution equa- 
tions for small perturbations. 




5F 



oc 



<S(A,S) 



M 



A 



(B2) 



where A represents deviations of the surface away from 
a pseudosphere, and S represents deviations away from 
an aster (^ = 0). Perturbations of the height field couple 
with perturbations of the angle field, so all four compo- 
nents of the two-by-two matrix of differential operators, 
M, are non-zero. The perturbations are functions of both 
radius and angle, and are generally not radially symmet- 
ric. To solve this, we write the perturbations in a Fourier 
basis. 



MT^ I A„(r) \ ^ime x^t 



(B3) 



where each two-vector (Am,Sm) is independent. Substi- 
tuting this solution into the evolution equation gives a 
separate set of coupled equations for each m-value. 

Neglecting bending rigidity, and choosing units of en- 
ergy such that a — \ and units of length such that 
K\ — 2a, we have for each value of m. 





dr-+r{l-r^)dl)) 


im{(2-K3)r-(2~(2~K3)r^)dr) ' 


l-r2 




im(if3r+(l-r^)(2-(2- 
r(l-r2)3/2 







(B4) 



For m = 0, the equations decouple. Since these are for S has such a solution, 
perturbations, we must find real-valued solutions that 

vanish at the boundaries r = and r — 1. The equation ^ gj^^ ( ) (B5) 




if Aq = —K^TT n for integer n. This is always negative. 
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The equation for A has the real-valued solution 

+r2SFi(&,6+^;26+|r2) , (B6) 

where 2F1 is the Gauss hypergeometric function 
and a = (-1 - Vl + 2Ao)/4 and 6 = (-1 + Vl + 2Ao)/4. 
Since the third argument exactly equals the sum of the 
first two, < r < 1 is the convergent domain for these 
functions. An ad hoc numerical study indicates that 
Ao —00 might extend this domain and allow the limit 
A(l) 0. These functions also diverge at r = 0, and 
again a large negative Aq appears to mediate this be- 
cause the function oscillates rapidly and might average 
to zero as r ^ 0. We lack an analytic treatment of this 
asymptotic regime, so we turn to a numerical method 
below. 

Considering < m and substituting u = r^, one sees 
that Eq. ( B4 1 consists of second-order ODEs with non- 



essential singularities at two points (the boundaries), so 
the equations can be transformed into hypergeometric 



differential equations |53j . By combining linearly inde- 
pendent solutions, one might construct real solutions that 
meet the boundary conditions for all values of m. Af- 
ter satisfying these constraints, one would obtain expres- 
sions for Am, which, when negative, indicate stable re- 
gions of parameter space. This approach is complicated 
even when rigidity is neglected. 

Instead of taking this approach, we have checked sta- 
bility numerically by discretizing the fields. We represent 
the deviations of the height and angle fields by a large 
column vector of field values at discrete steps in radius 
and polar angle. By representing the derivative operators 
as banded square matrices acting on this large vector, one 
obtains a matrix of numbers for any given set of parame- 
ter values. The largest non-zero eigenvalue of this matrix 
determines the stability of the shape. If the largest non- 
zero eigenvalue is negative, then that set of parameters 
suppresses perturbations and the shape remains stable. 

We have carried out such a numerical procedure. Gen- 
erally, the buckled aster is stable for 2k± < Ki < K3 and 
any < Kx- The analogous statement holds for buckled 
vortices. 
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